Subband structure of II- VI modulation-doped magnetic quantum wells 
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Here we investigate the spin-dependent subband structure of newly-developed Mn-based 
modulation-doped quantum wells. In the presence of an external magnetic field, the s-d ex- 
change coupling between carriers and localized d electrons of the Mn impurities gives rise to large 
spin splittings resulting in a magnetic-field dependent subband structure. Within the framework 
of the effective-mass approximation, we self-consistently calculate the subband structure at zero 
temperature using Density Functional Theory (DFT) with a Local Spin Density Approximation 
(LSD A). We present results for the magnetic-field dependence of the subband structure of shallow 
ZnSe/ZnCdMnSe modulation doped quantum wells. Our results show a significant contribution to 
the self-consistent potential due to the exchange-correlation term. These calculations are the first 
step in the study of a variety of interesting spin-dependent phenomena, e.g., spin-resolved transport 
and many-body effects in polarized two-dimensional electron gases. 
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Recent advances in nanoengineering of Mn-based semi- 
conductor heterostructures open up an exciting area 
of research which combines both magnetism and quan- 
tum size effects. The successful achievement of high- 
doping carrier densities in Mn-based II- VI layered com- 
pounds [|l] provides a new and promising ground for 
research: spin-polarized two-dimensional magnetotrans- 
ort in (semi)magnetic wells. Recently, Smorchkovaet al. 
reported Hall-effect measurements in spin-polarized 
two-dimensional electron gases in novel n-doped digital- 
magnetic heterostructures fi ||] . 

Here we investigate the magnetic-field dependence of 
the spin-dependent electronic structure of a modulation- 
doped digital magnetic quantum well (QW) || [|, q, 0, ||] 
using Density Functional Theory (DFT) with a Local 
Spin Density Approximation (LSDA). We consider a Dig- 
ital Magnetic Heterostructure (DMH) || in which mag- 
netic Mn atoms are digitally incorporated in a planar 
configuration inside the well. In the presence of an 
external magnetic field, the s-d exchange coupling be- 
tween carriers and localized d electrons of the Mn im- 
purities yields large spin splittings that exceed the inter- 
Landau level (LL) spacing ftw c ||. This gives rise to 
a magnetic-field dependent subband structure and a net 
spin-polarization of the electron gas even at low magnetic 
fields. 

Physical system. We consider the modulation-doped 
semi- magnetic quantum well shown in Fig. 1. The 
QW structure consists of ZnSe layers surrounding a 35- 
monolayer (ML) ZnCdSe/MnSe heterostructure with six 
equispaced digitally-incorporated MnSe monolayers, de- 
noted by 6 x i ML after Crooker et al. [|l0). This no- 
tation means there are 6 Mn planes with nominal pla- 
nar concentration x p — 1/2. Adjacent to each ZnSe 
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FIG. 1: Modulation-doped semimagnetic quantum well: the 
structure consists of ZnSe layers sandwiching a ZnCdSe layer 
containing digitally-incorporated MnSe monolayers. Adjacent 
to the ZnSe spacer there are n-doped regions (electron reser- 
voirs). Two depleted regions of length Ld arise within the 
n-doped layer as the doping electrons move into the well. 



spacer, there are n-doped regions (impurity concentra- 
tion Ni — 1.2 x 10 17 cm~ 3 ) that serve as electron reser- 
voirs. Electrons from the n-doped regions fill up the 
confined energy levels of the QW thus forming a two- 
dimensional electron gas (2DEG) and leaving behind de- 
pleted regions of length Lrj in the reservoirs. As the 
magnetic field is changed, we expect modifications in the 
self-consistent potential and in the corresponding band- 
structure due to i) magnetic-field dependent contribu- 
tions to the confining potential and ii) rearrangement of 
electrons in the confined levels because of the magnetic- 
field dependent Landau- level degeneracy eB/h. 

DFT/LSDA approach. Within Density Functional 
Theory, we solve the self-consistent Kohn-Sham equa- 
tions using a Local Spin Density Approximation scheme 
to obtain the magnetic field dependent subband struc- 
ture, modulation-doped confining potential, and carrier 
concentration. We consider the system at zero tempera- 
ture and perform the calculations within the framework 
of the effective-mass approximation. The dominant con- 
tribution to the magnetization comes from the localized 
moments of the Mn ions. These are treated as a param- 
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agnetic collection of ions by using an effective Brillouin 
function. However, we also consider the effects of the spin 
density on the exchange-correlation (XC) term. This is 
done within the interpolation scheme of Barth and Hcdin 
jnj as parametrized by Gunnarsson and Lundqvist 
Since the energy scale in our system (tens of meV) is 
the same as that of the calculated XC contribution to 
the confining potential (up to 18 meV), we expect XC 
to significantly affect the self-consistent subband struc- 
ture. Thus, to assess the role of exchange and correlation 
in the electronic structure in these novel geometries, we 
contrast our DFT/LSDA results with similar ones within 
the Hartree approximation Jl3| ]. 

Results and discussions. We present results for shal- 
low ZnSe/ZnCdMnSc modulation-doped QWs. By shal- 
low we mean the confining potential (25 meV) is small 
enough so that (<~) 10.5 nm wide wells have only a few 
confined subbands. For these geometries, n-doped re- 
gions with concentrations near 10 17 cm -3 provide carri- 
ers that fill up all the QW levels thus serving as a reser- 
voir of electrons ("reservoir hypothesis" 10). We obtain 
the magnetic-field dependence of the subband structure 
and our results show a significant contribution to the 
self-consistent potential due to the XC term. We have 
calculated the self-consistent potential profile, electronic 
structure and the two-dimensional electronic density for 
magnetic fields ranging from zero to ten tesla. 
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FIG. 2: Calculated modulation-doped potential profile, sub- 
band energy ei and wavefunction squared for spin-down elec- 
trons in the structure shown in Fig. 1 within the Hartree 
(dashed lines) and the LSDA (solid lines) approximations. 
The Fermi energy ef is set by the chemical potential of the 
reservoir, arbitrarily set to 25 meV. The main effect of the 
XC contribution in the LSDA calculation is to lower the con- 
fining potential and the subband energy Ei for the spin-down 
electrons when compared to the Hartree case. 



Figure 2 shows the potential profile, subband energy 
and wavefunction squared for spin-down electrons at B = 
1 T. The potential profile consists of the structural con- 
fining potential, the magnetic field and spin-dependent 
s-d exchange energy within the magnetic monolayers, the 



Hartree contribution, the Zeeman energy, and the spin- 
dependent XC energy. The resulting self-consistent po- 
tential has the usual modulation-doped profile that arises 
from the addition of the Hartree energy of the confined 
electrons to the confining potential. We contrast the 
LSDA calculation to the Hartree one. First we note that, 
despite the self-consistent Hartree result being fully con- 
verged, the LSDA energy profile is not flat at both ends 
as one would expect. The reason for this behavior may 
be that the LSDA parametrization of Gunnarsson and 
Lundqvist Q is not appropriate to describe XC at very 
low electronic densities such as those within the wave- 
function tails [their results are for Wigner-Seitz radius r s 
(in units of effective Bohr radius) within the range 1 — 9, 
while our data has r s from 20 to 300 at the tails]. Nev- 
ertheless, as most of the wave function is confined inside 
the well region, we do not expect these non-flat regions 
to drastically affect the electronic structure of the sys- 
tem. In fact, we have tested this hypothesis by abruptly 
setting the XC contribution to zero for r s greater than 9 
and we have not found any substantial deviation in our 
results. 




FIG. 3: Fan diagrams for the first three spin-down (tri- 
angle down) and spin-up (triangle up) Landau levels for 
both Hartree (a) and LSDA (b) calculations. Note that the 
exchange-correlation contribution enhances the spin-splitting 
(b) and makes the 2DEG fully polarized at lower magnetic 
field (~ 1 T) than in the Hartree case (~ 3 T). The Fermi level 
£/ (dashed line) remains constant at 25 meV (electron reser- 
voir) while the Landau levels oscillate as the spin-dependent 
potential profile changes with the magnetic field. 



We can clearly see from Fig. 2 that the XC contribu- 
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tion significantly lowers the energy E\ for the spin-down 
electrons. This effect is larger for the spin-down elec- 
trons than for spin-up electrons because the electron gas 
is spin-down polarized. Pauli's exclusion principle allows 
(does not allow) a spin up (down) electron and the major- 
ity spin down electrons to overlap thus raising (lowering) 
the electrostatic energy. 

In Fig. 3 we plot the first three spin-resolved Landau- 
level fan diagrams for both the Hartree (a) and the LSDA 
(b) cases. Note that the XC enhances the spin splitting 
because it decreases (increases) the spin down (up) en- 
ergy in the case of a spin-down polarized gas. The second 
point is that the LLs are pinned to the Fermi level e p for 
distinct ranges of magnetic fields (see plateaus at ef) [ fy§ - 
In the present model, this follows from the "reservoir hy- 
pothesis" which allows for two competing mechanisms to 
act: i) as the magnetic field increases, the LLs energies 
also increase and cross the Fermi energy, thus expelling 
electrons out of the quantum well; ii) on the other hand, 
less electrons means less Coulomb repulsion and the en- 
ergy tends to lower. An equilibrium between these two 
mechanisms is responsible for the plateaus. Note that 
the pinning of the Landau level in the present model is 
different from the standard one used to describe the In- 
teger Quantum Hall Effect (IQHE) [|J|. There, localized 
states in broadened Landau levels are responsible for the 
Fermi-level pinning Jl7t . 

One experimental relevant quantity affected by the 



electronic structure is the two-dimensional electron den- 
sity «2D- At zero magnetic field we find that n 2 D is 
2.175 x 10 n cm _2 and 1.444 x 10 n cm~ 2 , respectively, 
for LSDA and Hartree calculations. This considerable 
difference comes about because XC lowers the subband 
energies as compared to the Hartree case thus leading to 
an effective increase in n^D- 

We have studied the self-consistent electronic structure 
of a modulation-doped semimagnetic quantum well in- 
cluding exchange and correlation of electrons within the 
DFT/LSDA scheme. We found that for shallow geome- 
tries the exchange-correlation contribution is significant 
to the electronic structure as it i) lowers the subband 
energy of the majority spin-down electrons and ii) en- 
hances the spin-splitting thus making the electron gas 
fully spin-polarized at lower magnetic fields (~ I T) than 
in the Hartree case (~ 3T). Our results are preliminary 
but do point to the relevance of the XC in shallow ZnSe 
semimagnetic QWs, in contrast to GaAs systems. The 
calculations reported here are a first step in the study of 
a variety of interesting spin-dependent phenomena, e.g., 
spin-resolved transport and many-body effects in spin- 
polarized two-dimensional electron gases. 
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